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ABSTRACT 


The balance of wave action spectral density for a fixed wave-number is 
expressed in terms of a new dimensionless function, the degree of saturation, b, and 
is applied to an analysis of the variations of this quantity (and local spectral 
level) at wave-numbers larje compared to that of the spectral peak, that are produced 
by variations in the ocean surface currents in the presence of wind input and wave 
breaking. Particular care is taken to provide physically based representations of 
wind input and loss ty wave breaking and a relatively convenient equation (4.2) is 
derived that specifies the distribution of the degree of saturation in a current 
field, relative to its ambient (undisturbed) background in the absence of currents. 
The magnitude of the variations in b depends on two parameters, U 0 /c, where U 0 is th» 
velocity scale of the current and c the phase speed of the surface waves at the 
(fixed) wave-number considered or sampled by SAR, and S = (L/A) (u*/c)2, where L is 
the length scale of the current distribution, A the wavelength of the surface waves 
and u* the friction vr’ocity of the wind. When S is large 'of order 10 or more) the 
distribution of b is i r ensitive to currents for which U 0 /c ~1, but when S is of 
order unity or less, r gnificant variations in b are produced. A convergence zone is 
associated with a maximum in b relative to its ambient levels of 
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where m — 0.04 and n ~ 3. This would appear as a bright line in the SAR 
imagery. In general, local maxima in b (and the return SAR signal) should be 
observed if the local current strain rate scale 
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A local divergence or upwelling reduces the relative degree of saturation; when 
S is small the reduction is by the factor (1 + 2U/c) _ 9/c and continues until the 
waves grow back to the equilibrium level under the Influence c* the wind, A 
divergence line ould be Imaged as a line across which che return decreases from 
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the ambient level upwind, to a low ' level downwind, gradually recovering to the 
ambient. 

1 . Introduction 

The variety of instruments now available for active radio probing of the ocean 
surface makes possible the observation over a wide area of oceanic features ranging 
from capillary and microscale breaking waves to long ocean swells. Variations in the 
return signal with position on the ocean surface gives indications of even large scale 
features such as mesoscale ocean eddies and the distribution of surface wind stress. 

The scales of the surface features responsible for the return signal depends upon the 
wavelength of th incident radiation. When the angle of incidence on the water surface 
is not close to normal, and the wavelength of the radiation is large compared with the 
scales of occasional local occurrences such as breaking waves, the return signal seems, 
as a result of careful measurements such as those of Keller and Wright (1975) to be 
predominantly the result of first order Bragg scattering from freely travelling surface 
waves. The wave-number k of the surface wave sensed is 

k * 1*1 

where k^ is the wave-nuniber of the incident radiation and Q the angle of incidence; 
the depression angle is 1/2 IT- 6 . A real or synthetic aperture radar (SAR) at a fre- 
quency of 1.5 GHz responds to a surface wavelength of the order of 10cm at a depression 
angle of 30° or so, and a high resolution image such as that from SAR gives, in the 
pattern of intensity variations in the return signal, the pattern of local energy den- 
sity of surface disturbances at this scale. 

Some of these patterns are extremely intriguing. Beal et al. (1981) show a 
number of examples obtained from SEASAT, one noteworthy one giving surface expressions 
at this scale of tidal flow over the Davis Shoals, south and east of Nantucket Island. 
Another, off the mouth of the Chesapeake Bay suggests a complex pattern of eddies (or 
perhaps water masses); at high resolution the individual swell crests can be discerned. 
Fascinating as these images are, it is far from clear how they should be properly 
interpreted or what quantitave information they may contain. If we accept the premise 
that the images are r ived primarily from freely travelling short gravity waves, 
modified by longer waves and currents, then it should be possible to infer at least 
some properties of the current field, in particular, from the imagery. The short gra- 
vity waves sensed by the SEASAT SAR do have a significant ‘lifetime 1 of propagation, 
unlike capillary waves which are transient and fugitive, and their interaction with 
longer waves and currents may significantly modify their distribution over the ocean 
surface. The object of this paper is to determine the extent to which this can occur 
and to seek relations by which quantitative Information might be obtained from this 
imagery. 
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2. Action input and dissipation processes at short gravity wave scales . 

Since we will be concerned with interactions between the short gravity waves 
and surface currents, the wave dynamics are specified most conveniently by the balance 
of action spectral density N(k), which is related to the spectrum of surface displace- 
ment \y(k) by ~ 


N(jf) = ( ffrj '¥()£) = ^A) Vv 


( 2 . 1 ) 


cr 

where is the intrinsic frequency and*\^(k) is normalized such that 

j ■'$''(•*) ot -• S* , 

the mean square surface displacement. The spectral energy density of the wave field is 
|og\p(k), where p is the water density. 

Following energy paths, the balance of wave action spectral density is given by 

ft N(k) * J* * (C*t/)-VN -- - V k T(U) + S w - 0 , (2.2) 

(see, for example, Phillips, 1980) where C = represents the group velocity and U 
the velocity of the surface current. The~ exchanges of action density among different 
wave components by wave-wave interactions are specified by the first term on the right, 
in which T ( k, ) is the flux of action spectral density in the wave-number plane. These 
interactions, in a gravity wave system, are conservative so that the interval over all 
k vanishes. The two remaining terms specify the spectral distribution of action input 
from the wind and the loss of action density (which at these scales is predominantly 
the result of sporadic wave breaking or possibly the formation of parasitic 
capillaries). Our present degree of understanding of the three processes represented 
on the right hand side of (2.2) decreases monotonical ly in the order in which they are 
written though, a;', we shall see, some useful progress can be made. 

The specification of the wave field is completed by use of the klnematical 
conservation equation 
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together with the condition that • ° and the dispersion relation giving the 

intrinsic frequency in terms of the wave-number magnitude: 


<r * ( 9 k ) Vv 


(2.4) 


If the underlying current is non-uniform, the wave-number k of an energy packet 
varies as it propagates across the surface - k = k(x,t). On the other hand, in remote 
seising by radar (such as the L-band SAR), the Bragg scattering condition selects an 
almost fixed surface wave-number k in (1.1), so that it is convenient in this context 
to re-write (2.2) to refer to a fixed, rather than a variable wave-number. From (2.3), 
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but since Vxjc = 0 and Cj = i<r /ikj. 
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and since N = N(k,x,t), 
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( 2 . 6 ) 


now for a fixed wave-number k. 

To progress beyond this point, we need some assessment of the magnitude and 
nature of the terms represented schematically on the right of (2.2) or (2.5). The 
spectral energy and action transfers have been Investigated extensively by Hasselmann 
and his co-workers, by West (1984) and by Fox (1976). There Is little doubt that reso- 
nant spectral energy transfers Influence significantly the components near the spectral 
peak and, as suggested by Kltalgorodskll (1984) over a range of higher frequencies up 
to those for which c/u*^- 5 at most. Just above the frequency of the spectral peak, 
the fractional rate at which action Is acquired or lost by wave Interactions per wave 
period In an active wave field Is of the order 10~4 (Phillips, 1980) and this can be 
taken as an upper limit to the rate appropriate to smaller scales (0.1 to lm) of 
Interest here. Even so, it Is, as we shall see, generally less by at least an order of 
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magnitude than the rate at which action is acquired by snort gravity waves from the 
wind, provided the wind is sufficient to produce any such components at all. 

The energy transfer from wind to waves has been the subject of many theoretical 
and experimental enquiries during the past twenty-five years which have, if nothing 
else, demonstrated the complexity and variety of processes involved. At these small 
scales with wind blowing, the wave components are certainly of finite amplitude, with 
intermittently high curvature at their crests producing bursts of parasitic capillaries 
and occasional breaking. Air flow separation accompanies breaking as Banner and 
Melville (1976) have shown and Banner (1984) demonstrates that an enhanced local energy 
and momentum flux to the waves occurs as a result. Both the short waves and the wind 
are modulated and deflected by the longer waves and swell so that for the purposes of 
providing a more concrete expression for S w in (2.2), the best guide is given by the 
results of careful measurements interpreted in the light of only very general theoreti- 
cal consi H erations. Plant (1982), suggests from an examination of such measurements 
that 


S w (k) = 0.04«r( u»/e.)'cei9. Wfc) 


= tM U*/C ) X NCk 1 , 4*3, 


(2.7) 


where u* is the friction velocity, c = (g/k)*/2, the phase speed of the component and 
the angle between k and the wind. Plant believes the numerical coefficient to be 
accurate to within~about fifty per cent. Mitsuyasu and Honda (1984) give a similar 
expression derived from their experiments, with a numerical coefficient of 0.07 but 
with no directional factor included. Phillips (1980) estimated a coefficient of 0.05 
and Gent and Taylor's (1976) numerical calculations give approximately 0.07 — the 
relatively small scatter of values obtained from independent calculations or sets of 
measurements gives some confidence that this is a reasonably accurate representation 
over the range 1 > u*/c > 0.1 that includes most short wave cases of interest. The 
specific directional factor (cos 8) suggested by Plant is less certain. Nevertheless, 
the base of data he used is both extensive and carefully evaluated, so that (2.7) can 
be adopted with reasonable confidence as a useful semi -empirical representation. Since 
we are principally concerned with short wave components with wavelengths of 10cm or 
more and with wind speeds of order 10 m/s or less, the pertinent values of u*/c do 
occupy the range specified above and the fractional rate at which action is acquired 
from the wind per wave period. 
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ranges from 0.3 to 3 x 10-3 and is at least an order of magnitude larger than the net 
rate associated with wave-wave Interactions. 
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If the short wave components of the field have reached a state of statistical 
equilibrium in a region undisturbed by currents, then dN(k)/dt = 0 and since for these 
components V^.T is negligible, the densities of short wave breaking or parasitic 
capillaries must adjust themselves so that the dissipation D balances the wind input 
rate. Most previous observational studies on wave breaking such as those of and Toba 
and Kunishi (1970) have sought direct rela tionships between the density of whitecap 
coverage and wind stress. These implicitly assume such a statistical balance; our pri- 
mary concern here is the response of the short wave structure when the balance is 
disturbed. Whether or not a particular wave crest breaks depends, in fact, on its con- 
figuration and time development and not directly on the wind peed, since (for example) 
breaking can be induced by a local adverse current, independent of the wind sp.ed, and 
suppressed when the waves overtake a locally favourable current. Accordingly, D must 
be regarded intrinsically as a functional of the local spectral density, not the wind 
speed; its dependence on wind speed arises only insofar as this influences the spectral 
level or possibly, in the case of very short waves, as it influences the microscale 
breaking criterion through the surface drift (Phillips and Banner, 1974). 

The spectral level is conveniently specified by the dimensionless function 


B a H(k) c (2.8) 


which, in any particular spectnm, can be regarded as a function of k/k 0 , where k 0 is 
the wave-number of the spectral peak. This function can be called 'the degree of 
saturation'. The simple idea of a saturation range at high wave-numbers of the gravity 
wave (Phillips, 1958) assumed a hard upper limit B s , independent of k/k 0 , for B (the 
saturation range 'constant'), any further increase in the spectral level being prohi- 
bited by immediate wave breaking. Within the present conceptual framework, this 
corresponds to a dissipation function 0 that is negligibly small when B<B S but discon- 
tinuous (and in fact undefined) as B*B S . It is more realistic to assume that the 
action dissipation rate D is a continuous function of the degree of saturation B that 
is zero when B = 0, remains very small while B is significantly less than the old 
saturation range 'constant' when breaking is rare, and then increases rapidly when B 
approaches B s and short wave breaking becomes denser in space and more frequent in 
time. Its form is then as indicated schematically in figure 1. 

Also shown in figure 1 is the wind source term (2.7) which can now be written 
as 


) * Ooivco»6> (’u,/c.) v ^k’ M B(k/k # ) (2.9) 
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from ( 2 . 8 ). Fo r ' a particular wave-number th-‘s is linear in B to the accuracy of the 
expression ( 2 . 7 ) and the slope increases quadractical ly with u*. The short waves are 
in sat i st i cal equilibrium at the degree of saturation B defined by the intersection of 
these two curves. Note that, according to this description, the degree of saturation 
under equilibrium conditions may increase somewhat with wind speed, but the generally 
good observational support for the original idea of a 'hard' saturation range indicates 
that the variation wich wind speed must be weak and the curve of D as a function of B 
must increase very rapidly indeed. 

Before our specification can be completed for the spectral density of action 
dissipation, there is one further conceptual difficulty that must be faced. Does the 
function D for a particular wave-number depend on the degree of saturation at that 
wave-number only or on its value over a range of surrounding (or perhaps distant)' 
wave-numbers? An individual breaking event is localized in space within a distance 
probably of order k"l; in Fourier space its influence is distributed over a range of 
wave-numbers of order k. Also, the statistics of wave breaking at a certain scale 
(defined by the phase speed of the breaking crest) cannot be expected in general to 
depend on the spectral density or the degree of saturation at that scale alone, but 
over a range of scales. Consequently the statistical characteristics of wave breaking 
and the dynamical consequences of breaking on the wave field both involve a possibly 
wide spectral range. If, however, the degree of saturation B is almost constant over 
this range, a condition that might be anticipated under steady wind conditions, then on 
similarity grounds one might assert that 

D = gk-4 f (B) (2.10) 

In situations in which a much longer swell is present, the degree of saturation of the 
short waves will vary with respect to phase of the swell, being greater near the crests 
(with local breaking) and much smaller in the troughs. Although an expression such as 
(2.10) may hold locally, the overall mean rate of action di.slpatlon will depend not 
only on the mean degree of saturation but also on the swell slope, and this further 
parameter should be Included in (2.10). 

The function f(B) Is Illustrat'd only schematically In figure 1. Although Its 
detailed form Is not known at this stage. It may be possible to determine some of its 
characteristics more explicitly rrom laboratory or field observations. If the short 
waves are in a statistical balance between wind and dissipation. In the absence of much 
longer waves or swell, the degre.* of saturation attained Is defined by the balance of 
(2.9) and (2.10): 
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The variation of the degree of saturation with wind stress (normalized by phase 
velocity) is then 


*8 



( 2 . 12 ) 


According to the simplest ideas of a saturation spectrum independent of wind speed, the 
right hand side of (2.11) would be zero since the slope of f(B) is infinite at satura- 
tion. On the other hand, if a variation with wind speed of B as specified by (2.12) 
could be established reliably, then this would define the slope of the curve f(B) nea~ 
eqi.il ibrium. 

Another constraint is provided by the observation that the local action density 
of these short waves on the ocean is stable to perturbations on a scale large compared 
with the wavelength, in the sense that if locally the action density is high, the 
increased spectral density of breaking reduces it, and if it is low, the decreased 
spectral density of breaking allows the wind input to restore the equilibrium. The 
larger, dominant wave components at large fetch and duration may show a 'groupiness' as 
a result of the instabilities associated with non-linear wave interactions, but as we 
have seen, the time scales of these instabilities at higher wave-numbers are negligibly 
slow compared with those of the wind input and wave breaking. If under a constant wind 
field, the degree of saturation B is perturbed by an amount B', then with the neglect 
of the resonant interactions, we have for the rate of change of the perturbation in 
action spectral density N‘, 

— ' * c N (fc + a\) - DCS^a') 


Since N' = gk” 4 «‘“ 1 B ’ and in the equilibrium state S W (B) = D(B), this can be expressed, 
using (2.9) and (2.10) as 


* Lf<*) ~ , (2.13) 


where m » 0.04 cos 6 . 


For small perturbations B’ 
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and the condition for statistical stability is 




(2.14) 


which, in terms of figure 1, implies that the slope of the curve D as a function of B 
must be greater than that of S w at the point of intersection (as, indeed, we had 
anticipated). As u*/c increases, the equilibrium point rises and the slope df/dB at 
the equilibrium point increases even more rapidly. If the curve f(B) can be expressed 
reasonably accurately as a power law, the relaxation rate 


ir 




(2.15) 


can be expressed simply and explicitly 


in terms of u*/c. 


For, suppose that 
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f (B) = aB n 

where n is certai.ily areater than one, then at equilibrium 

r»(u+/ c.; v 8 - 

so that /Vaj * (u ,/ c.;*" 

and ir » -Of-Um/c)* o* 


(2.16) 


(2.17) 


Even though the numerical value of n Is not known (one might guess something in the 
range 3 to 5) this is a useful esult. It can be written alternatively as 




y * 
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(2.18) 


showing that, for a given wind stress, the relaxation rate towards equilibrium 
Increases as the cube of the frequency for these short wave components. 
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if, however, the perturbation in B is not small, the expected strong non- 
linearity in f(B) when B is small indicates that the recovery will be different for a 
local increase in degree of saturation than it will be for o decrease. If the decrease 
in B is such that f(B) is then negligible and breaking virtually ceases, the recovery 
will be initially at the undamped growth rate, (proportional to the slope of the line 
in figure 1 marked S w ) which may be much slower than * (proportional to the difference 
in the two slopes at the intersection.) A large local increase in B is then expected 
to disappear more rapidly by much enhanced breaking than a large local decrease in 
which breaking ceases. If, at equilibrium between breaking and wind input so that S w = 
D, the degree of saturation is B 0 , then from (2.9) and (2.10) 


m(u*/c)2 gk' 4 B 0 = gk" 4 aB 0 n 

when the representation (2.16) is used. Consequently 

a r m(u*/c)2 Bo 1 " 0 

and the balance of the last two terms on the right hand side of (2.16) can be written 
in terms of this equilibrium level as 


gk- 4 . m(u*/c)2 (1 - (B/B 0 )n-I)B 


(2.19) 


3. The influence of a background sea state on mean short wave properties at given k . 

The short ocean wave components that are generally responsible for radar back- 
scattering do not, of course, exist in Isolation. Longer waves or swc 1 produce modu- 
lations In the Intensity of the back-scattered return, measured by "right et al (1980) 
w. h a spatial scale of variation that Images the long wave field. Be 1 (1981) has 
made use of these modulations to study the evolution and propagation of swell across 
the Atlantic continental shelf of the United States. Frequently, the orbital veloci- 
ties associated with these longer waves are larger than the propagation speeds of the 
short components sensed and ones Intuitive feeling might be that an Irregular or random 
field of longer waves would might so scramble the propagation characteristics of the 
shorter waves as to obscure any likelihood interpreting the large scale patterns In 
their Intensity In terms of currents. However, this particular concern does not seem 
to be justified. 

Consider a surface current field u(x,t) whose scales of spatial and temporal 
variations L and T are large compared with the wavelength and period of long waves or 
swell, the scales over which the tangential components of their orbital velocity 


442 


ORIGINAL Page ra 

* p °o« o sag 


varies. Local mean values of the wave field, indicated by an overbar, are defined over 
scales that are small compared with those of the current field but large compared with 
those of the swell. Fluctuations or modulations produced by the swell are represented 
by the difference between the actual values and the mean values defined in this way. 
Thus, the action density is represented as N(x,t) + N 1 , ;e short wave group velocity a 
constant k_ (which varies because of the variations in th«- local g at the surface as 
the short waves ride over the swell) by C + C‘ , and so forth. If the short waves are 
everywhere close to saturation under the balance of wind input and wave breaking then 
(2.5) can be written as 
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(3.1) 


where resonant wave-wave interactions are neglected, N 0 is the equilibrium value of N 
in the absence of long waves or currents and r the relaxation rate of the previous 
section. Substitution of local mean and fluctuating quantities into (3.1) and 
averaging gives 
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The left hand side of this equation has the same form as the original equation (2.5) 
so that the influence of the swell on the mean action spectral density Is represented 
by the two new terms on the right. Now, to the first order, the wave-induced modula- 
tions N' are proportional to N (Wright et al 1980) and to the root-mean square slope 
of the long wave field, s are the local gradients with respect to x. and k^ , so that 
the covariance terms are of ordered* , whereas the swell frequency. Their only 
effect, then, is to modify slightly (by order e* ) the equilibrium mean spectral levels 
from that which would obtain in the absence of swell. 

The fact that the swell can therefore be Ignored In a consideration of the 
local mean back-scattering from short propagating waves Is at first sight counter- 
intuitive, but the essence of the matter Is that the modulation patterns at a fixed 
wa ve-number are periodic and non-cumulatlve provided the swell slope Is sufficiently 
small that the dissipation response can be linearized. 
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4. Variations In N produced by surface currents . 

The dynamical balance (2.5) can be written In term;* of the degree of satura- 
tion B by substituting n * lc' v *. 0(k) and the expression (2.18). 
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If the magnitude of the current variations Is represented by U 0 and the length scale 
over which they occur by L, then 


U 1 « Uo fl(x/L) - U 0 f(l), say 


Also, If b(T) * B/B 0 , the local degree of saturation relative to Its equilibrium 
level, then~under steady conditions (4.1) can be written In dimensionless form as 



q U, <Cc«C i 
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where c represents the phase velocity and S the (fixed) wavelength of the compel 
considered. 

It Is Interesting to note that this equation nvolves only two basic , 'pe- 
ters, U 0 /c expressing the strength of the current field and (L/A) {cj c)2 repre- 
senting the combined Influence of the wind, the scale of the current field and. 
Implicitly, breaking. This second parameter can also be written as 


fj 




0«r as 




444 


ORIGINAL PAGE fS[ 
OF POOR QUALITY 


where C is the group velocity. The latter form provides a useful physical interpreta- 
tion of this parameter as proportional to the ratio of the time taken for an energy 
packet to move over one scale distance L to the characteristic growth time of the 
waves under the influence of the wind alone. If this para meter is large, re- 
adjustment of the waves is rapid in the time taken to tra verse the current variation, 
so that the resptnse in the relative degree of saturation is small. Small values of 
the parameter indicate situations where the re-adjustment time is long, so that one 
might expect a significant response. 

For the sake of definiteness, suppose that the 1-direction in (4.2) is chosen 
as the direction of the wave-number being sampled. The middle term on the left-hand 
side then reduces to 


4 u. it 

2 . c ar, 


expressing the influence of convergence or divergence of the current field in this 
direction. The third term involves the variation with respect to magnitude and direc- 
tion of the degree of saturation and must be considered when a current, varying In x, 
is in the transverse direction so that energy paths are deflected. Little is known 
(although much is speculated) about the angular distribution of the wave spectrum cr 
b, so that in order to obtain a sense of the magnitude to be expected of the 
variations in the local degree of saturation relative to the background, we will 
neglect the variations uf b with k ano consider convergences or divergences only with 
* U(xi). In the flov over r’atively shallow topography such as the Nantucket 
Shoals, both local convergences (iU/ix < 0) and divergences (2u/i x > 0) can be 
expected; under open sea conditions, divergences would be expected in regions of local 
upwelling whle local convergences may be found along fronts, with one water mass 
overriding another, slightly denser. 


With f(t) i ( f(r) t o, o ) , say, equation (4.2) reduces to 
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where S Is the parameter (L/ > )(u # /c) 2 . The value of the Index n, although greater 
than 1, Is not ki.own at this time; for the purposes of these calculations n was taken 
as 3. A smaller value gives a more gradual onset of energy loss by breaking as b 
Increases and so a larger response to convergences; a larger value reduces the 
response. This is shc^n explicitly below. The wind is supposed to act In the x 
direction also so that m - 0.04 and a current distribution of the form 


f(l) » i f » ♦ f} 
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was chosen. The equation (4.3) was integrated numerically with b = 1 at S = -2 for 
various values of U 0 /c and S, using an Apple II computer. When U 0 /c < 0.5, equation 
(4.3) has a mathematical singularity at which the coefficient of the first turn 
vanishes. In those cases, the integrations were performed following the direction of 
the wave energy flow, from S - -2 to within o.ie bteo o* the singularity and then from 
S= + 6 bach to the singularity iince in this part of the physical domain, the backwards 
energy convection of the current overcomes the forwards propagation. If the waves on 
the far side of the convergence zone have had sufficient fetch and duration to achieve 
equilibrium, then b = 1 at J = 6 also and it is found, somewhat surprisingly, that the 
two branches of the solution form a continuous curve across the singularity at which 
(1/2 + (U 0 /c) f ( S)) 3 0. On the other hand, if the fetch or duration on the far s^e 
is limited, then b on the far side may be less than unity for these waves and a real, 
physical singularity will occur a :his point. It seems that only ir this casr will 
an abrupt change in the return signal coincide with this singularity. 

Seme results of these calculations are illustrated in Figure 2. Vhe n the 
parameter S is of erder 10 or more, the responses of b to the current field are srnal’ 
and 11mite J essentially to the region of current shear. Fir smaller values of S, the 
magnitude of the response Increases; in a convergence (U 0 /: < 0) the region ova which 
b is significantly greater than 1 remains limited to the onvergence region because of 
the relatively rapid reduction In wave energy by breaking. In a divergence, however, 
the degree of saturation at a given V. reduces as a result of two effects - the 
divergence spreads the action density out over a greater spatial Interval and (more 
important) the energy packet at the wave number observed after divergence was, before 
the divergence, at a larger wave-number or smaller scale, so that its action density 
was less. At :»lues of S less than aoout 4, the recovery distance of the wave field 
after suppression is significantly greater than L sc tnat the waves remain unsaturated 
for a cons1de~able distance. 

In a convergence zone, the maximum degree of saturation relative to the 
background occurs close to the point of maximum strain rate, when f'(J) » 1. From 
(4.3) this occurs when 

- J S (i- ) % 


so that the maximum contrast Is 

* I * +■ n ) (4.4) 


In the cases shown In Figure 2, n been taken as 3; since the coefficient 9/4«m » *o 
20, hoax varies approximately as (*i 0 /c) S-% for values of S < 4 U 0 /c or so. 8ut the 
value of the index n In the representation (2.15) cannot be regarded as known and. In 
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fact, the expression (4.4) may aiiow it to be determined by field measurements of the 
maximum contrast in the degree of saturation at a fixed k in convergence zones over a 
•'ange of values of S and U 0 /c, although if n is significantly greater than 1 (as 
expected), the dependence on n is rather weak. 

A general conclusion from these computations is that if freely propagating 
waves under the influence of wind and breaking are to experience a significant 
variation in the degre» .1 saturation at a fixed wavelength A , then the second term 
in the bracke. of (4.4, must be greater than unity, say, or that 


U * 
c S 


o os 


(4.5) 


Expresses alternatively in terms of the strain rate in a convergence, 
response requi res 


Uo 

u 


£ 0-12. 


ui 




a significant 


(4.6) 


For example, if u* = 20cm/s and A = 16 cm, then U 0 /L must be greater than about 2 x 
10“2 sec-1, which is a very large oceanic strain rate outside frontal regions. Light 
winds clearly favour a more sensitive response. 


It should be noted, however, that for a given convergence field and wind 
stress, (4.4) can be written as 


km* « 




U u? i > 


(4.7) 


the contrast in degree of saturation increases with wavelength, the smaller value of 
U 0 /c being more than compensated by the increase in S-l. Increased harmonic content 
of these longer waves and denser breaking patches may give enhanced Bragg and specular 
scattering that may possibly dominate the direct return from freely travelling waves 
at the wave-number sensed by Bragg scattering in the background. 


In a region of surface divergence (and upwelling), when S is small the reco- 
very of the wave field by the wind is slow over the scale of the divergence, and the 
reduction In b can be estimated by neglecting the wind input (and wave breaking) over 
this interval. With U 0 f(f) = U(T), equation (4.3) reduces to 


a - % r> ) f, 


5 ^ 

i ai 


( % v ) 


<r ~ ^ 
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which integrates exactly to 




t>(fj r (it xu(f)/c ) > (4.8) 

since b = 1 when U = 0. This reduction (U(X) > 0) at a fixed wave-number is, of 
course, greater than for an energy packet or wave train at a fixed apparent frequency, 
as calculated by Longuet-Higgins and Stewart (1960) since the waves elongated 
by the divergence to the wave-number k had a larger wave-number and smaller spectral 
density than the waves at wave-number k before the straining. 

The SEASAT synthetic aperture radar imagery, such as that reproduced in the 
book by Beal et al (1981; see especially pp. 22 and 96), offers a number of examples 
that seem, in the light of these results, to be interpretable as lines of local con- 
vergence or divergence. There is unfortunately little independent information on the 
current field or the local wind stress so that attempts at a quantitative comparison 
must av.ait imagery for which more complete documentation is available. Hovever, these 
results do indicate the potential in synthetic aperture radar imagery for measuring 
rates of strain in oceanic fronts and rates of upwelling over a substantial spatial 
domain. 
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Figure 1: Input from the wind, S , and dissipation D of action 

spectral density as functions of the degree of saturation 
B. The intersection specifies B Q , the high wave-number 
equilibrium value. 
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Figure 2: Relative degrees of saturation b * B/B q at a fixed 

wave-number as functions of position in a local convergence 
or divergence centered at x “ 0. In Figure 2a, 

(L/X)(u*/c) 2 - 12; in 2b, 4; in 2c, 1 and in Figure 2d, 

0.5. The three upper curves in each group represent the 
convergence cases with U Q /c ■ -1.2, -0.8 and -0.4; the 
three lower curves are for divergences with U Q /c - 0.4, 

0.8 and 1.2. The black dots represent positions where 



